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Nomenclature  (continued) 


<j>  angle  between  tangent  to  body  surface  and  xq  -  direction 

* 

p  density  (dimensional) 

f00 

0  momentum  deficit  area  =  — (1  -  — )rdy 

*  0  e  e 

x  shear  stress  in  boundary  layer  (dimensional) 

* 

wall  shear  stress  (dimensional) 


All  other  quantities  are  defined  in  the  text. 

All  quantities  in  the  text,  unless  otherwise  specified,  are  made 
dimensionless  as  follows: 

distances  with  respect  to  L 

i  .  * 

velocity  with  respect  to  U 

00 

"k  *2 

pressure  with  respect  to  p  u 
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Introduction 

This  report  addresses  the  problem  of  improving  the  detailed  predictions 
of  thick  axisymmetric  turbulent  boundary  layers  using  a  standard  finite 
difference  calculation  scheme  and  a  two-piece  algebraic  turbulence  model. 

The  typical  result  of  such  calculations  is  that  in  the  vicinity  of  the 
tail  the  mean  velocity  profiles  are  too  full  close  to  the  wall  which  is 
a  reflection  of  the  Reynolds  stress  being  much  too  high  across  the  boundary 
layer.  The  cause  of  this  erroneous  behavior  is  the  use  of  a  turbulence 
model  which  accounts  for  only  the  thin  boundary-layer  rate  of  strain 
proportional  to  3u/9y  and  neglects  the  extra  rates  of  strain  due  to  longi¬ 
tudinal  and  transverse  curvature.  The  experimental  work  of  Patel 
and  Lee  [1]*  has  shown  that  these  two  extra  rates  of  strain  are 
extremely  important  in  determining  the  turbulence  properties  in  a  thick 
axisymmetric  boundary  layer.  They  have  also  presented  a  method  for  the 
incorporation  of  these  extra  rates  of  strain  in  a  one  equation  turbulence 
model  and  made  calculations  with  their  model  which  show  excellent  agreement 
with  experiment.  In  this  report  the  extra  rate  of  strain  approach  of  Patel 
and  Lee  is  modified  for  application  to  a  standard  two-piece  algebraic  turbu¬ 
lence  model. 

Turbulence  Model 

The  Reynolds  stress  is  treated  by  a  simple  mean  field  closure  which 
in  a  thin  boundary  layer  reduces  to 

-  u*  v*  *  p*e*  (1) 


*Numbers  in  brackets  denof  ■»  Referenr  s  cited  at  the  end  of  the  report. 
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where  e*  is  the  eddy  viscosity.  A  zero-equation,  or  algebraic,  model  is 
«  used  where  the  eddy  viscosity  is  related  directly  to  the  mean  velocity. 

Such  a  model  requires  the  boundary  layer  to  be  divided  into  an  inner  (wall) 
region  and  an  outer  (wake)  region  with  a  different  equation  for  e  in  each 
region,  as  a  consequence  of  the  two  different  length  scales.  The  junction 
point  between  these  two  regions  is  taken,  in  the  usual  way,  to  be  the  value 
of  y  where  the  inner  and  the  outer  values  of  £  are  equal. 


With  the  aid  of  Eq .  (1),  the  shear  stress  in  a  boundary  layer,  which  is 
the  sum  of  molecular  and  turbulent  contributions,  can  be  expressed  as 


T*  =  (W*  +  0*£*)  ^  . 


(5) 
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Equation  (10)  is  a  generalization  of  the  form  discovered  by  Rao  which 

reduces  to  his  expression  when  <p  =  0  (when  the  body  is  a  cylinder  or  when 

4 

it  has  a  point  of  zero  slope) .  In  the  form  used  by  Cebeci  for  Y  ,  cos  <}> 
is  omitted.  Strictly  speaking,  Eq.  (10)  is  valid  only  for  small  pressure 
gradient  and  hence  4>  cannot  be  very  large  (say,  less  than  0.2). 

Since  Y+  formally  transforms  a  thick  axisymmetric  boundary  layer  into 
a  two-dimensional  one,  we  can  obtain  the  eddy  viscosity  for  the  thick 
axisymmetric  case  by  substituting  Y+  for  y+  in  the  two-dimensional  expression. 
Thus  for  a  thick  axisymmetric  boundary  layer  is 


2  +  + 
z.  ■  a* 

1  r  dy 

o 


(11) 


The  mixing  length  £  in  the  viscous  dominated  region  near  the  wall  as 
well  as  in  the  law-of-the-wall  region,  is  given  by 


£+  =  0.41Y+  D, 


(12) 


where  D  is  the  Van  Driest  damping  function  which  supresses  the  linear 

4  4 

dependence  of  £  on  Y  near  the  wall.  The  expression  for  D  is 

D  =  1.0  -  exp(-Y+/A+)  ,  (13) 

where  A+  is  an  effective  sublayer  thickness,  sometimes  referred  to  as  the 
Van  Driest  damping  parameter.  The  parameter  A+  determines  the  thickness 
of  the  sublayer  and  is  a  function  of  the  pressure  gradient  and  transpiration 
rate. 


Sublayer  Damping  Parameter 

The  expression  used  here  for  A+  is  the  one  deduced  from  experimental 
data  by  Crawford  and  Kays  [5],  as  follows: 


26.0 


r 

a  i  v  +  b 
i  w 


A 


:  I 

\l+cv  ) 


+  1.0 


(14) 
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Then  Eq.  (4)  becomes 


r*(y*  +  p*e*) 


9u*  - 
-5— r  =  r*T* 
9y*  o  w 


which  in  universal  law-of-the-wall  variables  is 


r  dy 

o 


=  1 


(6) 


The  thin  boundary-layer  counterpart  of  Eq .  (6)  is 


<1  +  0  *■_  - 1 

3y 


(7) 


Rao  [4]  was  the  first  to  show  experimentally,  on  a  cylinder  with  its  longi¬ 
tudinal  axis  aligned  with  the  flow,  that  the  law  of  the  wall  still  holds 
provided  the  normal  spatial  variable  is  redefined  appropriately.  Mathematic¬ 
ally  Rao's  discovery  means  that  Eq.  (6)  can  be  transformed  into  Eq .  (5)  by 
a  change  of  independent  variable.  If  Y+  denotes  the  new  axisymmetric  variable, 
then  we  require 


JL  =  JL  Ll  =  1 _ r+ 

9y+  9  y+  dy+  9y+  r  + 


and  hence, 


dy+  r+ 

dY+  "  r  + 
o 


(8) 


From  the  geometry  of  a  thick  axisymmetric  boundary  layer,  as  shown  in 

Fig.  1. 


r  =  rQ  +  y  cos  . 


(9) 


Equation  (9)  allows  Eq .  (8)  to  be  integrated  yielding  the  axisymmetric 
law-of-the-wall  variable 

_  +  _+ 


Y+  = 


cos  <p 


Jin 


(10) 
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(15) 


(16) 


(17) 


Equation  (14)  has  been  derived  for  boundary  layers  in  near  equilibrium 

■f  + 

where  P  and  v^  vary  at  most  slowly  along  the  surface. 

When  rapid  changes  occur  in  P+  or  v  ,  the  sublayer  does  not  adjust 

+ 

instantaneously  to  the  new  conditions.  Thus  A  will  lag  its  equilibrium 
+  4-4* 

value.  Since  A  is  a  function  of  P  and  v  ,  Crawford  and  Kays  have 

w 

introduced  lag  equations  for  these  parameters  to  allow  A  a  finite  time 
(or  distance)  to  adjust  to  a  new  equilibrium  state.  The  lag  equations  are 


of  the  form 


(18) 


where  denotes  either  P+  or  v+  and  the  subscript  e,  the  equilibrium  value. 

w 

Thus  $+  is  the  effective  value  of  P+  or  v+  to  be  used  in  Eq .  (14).  The 

w 

quantity  C  is  a  lag  constant  with  the  recommended  value  of  4000,  according 


to  Crawford  and  Kays. 
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Unfortunately,  for  a  body  of  revolution  Eq.  (18)  gives  erroneous 


behavior  near  the  tail.  What  happens  is  that  d$+/dx+  changes  sign  causing 
$  to  diverge  from  which  is  physically  unreasonable.  The  cause  of 
the  sign  reversal  is  the  use  of  x+  as  an  independent  variable.  By 


definition, 

x+  =  Re  u^  x  .  (19) 

As  the  tail  is  approached  u^.  begins  to  fall  rapidly  to  zero,  as  it  must, 

+  -f* 

so  that  x  reaches  a  maximum  and  decreases.  At  the  maximum  of  x  ,  the 

derivation  d$+/dx+  changes  sign.  This  erroneous  behavior  can  be  easily 

corrected  by  the  following  slight  modification  of  the  lag  equation: 


+  <t+  -  $+ 

d$  e 


-  (20) 


where  C(x)  = 


Re  uT(x) 


Equation  (20)  contains  the  same  characteristic  length  scale  as  Eq .  (18) 
but  because  x  is  used  as  the  independent  variable  no  sign  reversal  occurs. 

The  two  lag  equations,  (18)  and  (20),  have  been  tested  on  a  modified 
spheroid  at  a  Reynolds  number,  based  on  chord  length,  of  1.262  x  10^  with 
=  0,  corresponding  to  the  experiment  of  Patel,  Nakayama  and  Damian  [6]. 
Boundary-layer  calculations  were  made  using  a  Cebeci  and  Smith  type  of 
finite  difference  code  developed  at  ARL.  The  measured  pressure  distribution 


was  used  in  the  calculation.  For  this  body  x  has  a  maximum  which  occurs 

at  x  =  0.77.  The  variations  of  P  and  P  with  x  are  shown  in  Figs.  2  and  3 

e 

for  Eqs.  (18)  and  (20)  respectively.  The  choppy  character  of  P+  exhibited  in 

these  figures  is  caused  by  numerical  inaccuracy  in  computing  dp/dx.  Figure  2 

"f1 

shows  clearly  the  divergence  of  P  from  P  which  occurs  downstream  of  x 

e  max 
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when  Eq.  (18)  is  used.  When  Eq.  (20)  is  used,  as  shown  in  Fig.  3,  P 
lags  P+  in  a  plausible  manner.  We  would  expect  the  wall  friction  coefficient 
to  be  a  sensitive  measure  to  differences  in  A+  and,  hence,  in  p"*".  Figure  4 
reveals  that  the  calculated  values  of  c^  using  Eqs.  (18)  and  (20)  are 
indistinguishable  until  x  becomes  larger  than  about  0.86.  Beyond  x  =  0.874 
the  boundary-layer  calculation  using  Eq.  (18)  failed  to  converge.  As  can  be 
seen  in  Fig.  2,  at  x  =  0.874,  P+  becomes  negative  and  decreases  very  rapidly 
An  alternate  expression  for  a"*"  has  been  derived  by  Cebeci  [2]  as  an 
extension  of  Van  Driest's  one-dimensional  unsteady  analogy  and  is  given  by 

A+  =  &  ,  (22) 

N 

where  [1  ~  exp(11.8v  +) ]  +  exp(11.8v+)  .  (23) 

■+■  w  w 

V 

w 

+  ”4* 

Cebeci  does  not  use  lag  equations  for  P  or  v^. 

For  the  case  of  v+  =  0,  considered  herein,  the  variation  of  A+ 
w 

with  P+  according  to  the  empirical  expression  of  Crawford  and  Kays,  Eq .  (14), 
and  the  analytical  expression  of  Cebeci,  Eq.  (22),  is  shown  in  Fig.  5.  The 
two  expressions  differ  markedly  for  P+  t  0  with  the  differences  being 
largest  for  PH  <  0  (favorable  pressure  gradient).  Because  Eq .  (14)  is  a 
fit  of  experimental  data,  it  was  chosen  for  use  over  Eq.  (22). 


Outer  Eddy  Viscosity 

In  the  outer  portion  of  the  boundary  layer  (law  of  the  wake  region)  we 
have  the  choice  of  using  a  constant  mixing  length,  as  used  by  Crawford  and 
Kays,  or  a  constant  eddy  viscosity,  as  used  by  Cebeci.  The  choice  between 
the  two  in  the  present  instance  was  dictated  by  computational  reasons. 

The  matrix  associated  with  the  finite  difference  form  of  the  boundary- 
layer  equations  at  a  particular  x-wise  station  becomes  ill  conditioned  near 
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the  outer  edge  of  the  layer  for  the  constant  mixing  length  formulation. 

When  a  linear  solver  without  pivoting  is  used  the  resulting  velocity 
profile  exhibits  fluctuations  (numerical  noise)  which  make  the  location 
of  the  outer  edge  very  difficult  to  determine.  With  the  constant  eddy 
viscosity  formulation  these  fluctuations  are  absent  which  was  the  reason 
it  was  chosen.  Later  on  the  discovery  was  made  that  by  using  a  solver 
with  pivoting  (which  nearly  doubles  the  computation  time  per  profile)  the 
fluctuations  in  the  velocity  profile  for  the  constant  mixing  length 
formulation  were  effectively  eliminated. 

In  the  outer  region  of  the  boundary  layer  the  eddy  viscosity  is  taken 
to  be  constant  using  the  form  first  suggested  by  Clauser  and  later  modified 
by  Cebeci  for  low  Reynolds  number  effects  [2]. 

e  =  Re  a  u  <5*  ,  (24) 

o  e  k 

■jlf 

where  the  subscript  o  denotes  "outer"  eddy  viscosity,  and,  6^  =  kinematic 
displacement  thickness 


j»00 


and  a  is  a  dimensionless  constant  given  by 


where 


i  +  n 

a  =  a0  mr 

a  =  0.0168 
o 

n  =  0.55 

0 

n  -  n  [1  -  exp(-0.243  vrZ^  -  0.298Zj)] 


Re 


0k 


425 


(25) 


(26) 

(27) 

(28) 
(29) 


1 


(30) 


-16- 


30  September  1982 
GHHrcfb 


Some  formulations  also  include  an  intermittency  factor  to  account  for  the 
intermittent  character  of  turbulence  at  the  outer  edge  of  the  boundary  layer. 
In  the  present  formulation  this  factor  is  omitted  as  it  has  been  found  to 
have  negligible  effect  on  the  velocity  profile  and  other  boundary-layer 
parameters. 


Curvature  Effects 

The  so-called  thick  boundary- layer  formulation,  discussed  in  the  section 
on  the  inner  eddy  viscosity,  although  necessary  in  defining  the  correct 
boundary-layer  thickness,  gives  inadequate  results  for  the  details  of  the 
boundary  layer  as  the  tail  of  an  asymmetric  body  is  approached.  The 
predicted  Reynolds  stress  is  always  much  too  large  and  the  resulting  velocity 
profiles,  therefore,  much  too  full.  Many  investigators  have  noted  this 
deficiency  [1,  6,  7]. 

This  defect  in  axisymmetric  boundary-layer  calculations,  as  mentioned  in 
the  introduction,  may  be  traced  to  the  neglect  of  the  extra  rates  of  strain 
due  to  transverse  and  longitudinal  curvature  which  produce  important  extra 
terms  in  the  Reynolds  stress.  A  good  bibliography  of  the  experimental  work 
on  the  effects  of  curvature  on  turbulence  as  well  as  a  discussion  of  the 
dynamics  qf  the  problem  may  be  found  in  the  report  by  Patel  and  Lee  [1]. 
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Bradshaw  [8]  proposed  a  relation  to  account  for  the  effect  of  extra 
rate  of  strain  in  a  thin  shear  layer.  His  proposal  was  a  linear  modifica¬ 
tion  of  the  turbulent  length  scale,  valid  for  small  extra  rate  of  strain, 
as  follows: 


L  =  i+% 

£  3u 

°  3y 


(33) 


where  £q  is  the  length  scale  (mixing  length)  associated  with  the  rate  of 
strain  3u/9y,  £  is  the  length  scale  with  extra  rate  of  strain  e  included 
and  S  is  a  constant  of  order  10.  Bradshaw  also  proposed  that  the  equilibrium 
value  of  e  in  Eq.  (33)  be  replaced  by  an  effective  e  to  account  for  rapid 
changes  in  the  boundary-layer  history.  The  effective  e,  denoted  by  ee£^» 
is  then  governed  by  a  lag  equation, 


d_ 

dx 


(eeff)  L  (e  "  eeff) 


(34) 


e 

where  Le  is  a  lag  length  for  the  boundary  layer  response  to  changes  in  e. 
The  suggested  value  of  Lg  is  106. 

Bradshaw  recognized  that,  in  general,  e  would  be  a  function  of  y  as 
well  as  x.  To  get  around  this  difficulty,  since  Eq.  (34)  is  an  ordinary 
differential  equation  in  x,  he  suggests  removing  the  y-dependency  in  e  by 
using  e/u  in  Equation  (34)  in  cases  where  e  a  u  such  as  curvature,  lateral 
divergence  and  bulk  compression. 

Patel  and  Lee  [1]  have  incorporated  Bradshaw's  ideas  concerning 
curvature  effects  into  an  axisymmetric  finite  difference  boundary  layer 
calculation  scheme  with  a  one-equation  transport  model  for  the  turbulent 
shear  stress.  Comparisons  of  their  own  extensive  data  for  the  two  body 
shapes  with  different  curvature  histories  with  the  calculations  show 
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very  good  agreement  when  both  longitudinal  and  transverse  extra  rates  of 
strain  are  included.  The  extra  rates  of  strain  in  these  calculations  are 
.  not  small  so  that  Bradshaw's  linear  formula  breaks  down,  leading  to  a 
very  rapid  decrease  in  £  across  the  boundary  layer  and  thus  to  a  near 
total  destruction  of  the  Reynolds  stress.  To  correct  this  situation  they 
propose  the  following  nonlinear  formula: 


L 

£ 


a  e 


1  - 


ef  f 


o  — 


9u 

3y 


-1 


(35) 


which  reduces  to  Bradshaw's  linear  formula  when  a  e  ,,./u  is  small. 

eft  y 

Huang  et  al.  [7]  have  suggested  that  the  extra  rates  of  strain  in  a 
thick  axisymmetric  boundary  layer  can  be  accounted  for  by  the  following 
simple  formula: 


1/2 


L 

£ 


[(1  +  6/ rQ)i  -  11 

3.33  6/r 


(36) 


This  formula  is  based  on  the  idea  that  the  square  of  the  turbulent  length 
scale  is  related  to  the  annular  area  of  the  boundary  layer.  Huang's 
formula  gives  a  uniform  reduction  of  the  two-dimensional  mixing  length 
across  the  entire  boundary  layer  therefore  producing  too  large  a  reduction 
in  Reynolds  stress  in  the  wall  layer  region. 


Incorporation  of  Curvature  Effects  in  the  Turbulence  Model 

The  formula  of  Huang  as  well  as  a  number  of  variations  of  the  nonlinear 
formula  of  Patel  and  Lee  were  incorporated  in  the  present  two-piece  algebraic 
turbulence  model  and  applied  to  several  axisymmetric  bodies  with  different 
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curvature  histories.  The  best  formulation,  as  judged  by  agreement  between 
the  predicted  results  and  experiment,  was  found  to  be  the  following  modifi¬ 
cation  of  the  formula  of  Patel  and  Lee: 


j~  «  (1  -  f$)_1 


where: 


f 


A 


a 


3y 

and  e^  =  transverse  extra  rate  of  strain 

=  _ u _  1  dr 

et  1+ky  r  3x 

e^  =  longitudinal  extra  rate  of  strain, 
_  ku 

e£  1+ky  * 


(37) 


(38) 

(39) 


(AO) 


(AO) 


From  Eq.  (9)  and  the  definition  of  body  longitudinal  curvature  K,  we  find  that 


~  dr 
dr  _  _ o 

3x  dx 


-  <y  sin  <}>  , 


and  hence  Eq.  (39)  can  be  written  in  the  final  form: 


$  =  sin  *  -  r]  •  <42> 

3y 

The  manner  of  implementation  is  as  follows: 

1.  The  extra  rate-of-strain  effects  are  gradually  switched  on  beginning  at 
x  =  xq,  the  station  where  6  >  0.1,  corresponding  to  where  the  thin 
boundary- layer  formulation  begins  to  break  down. 

2.  The  exact  formulas  for  e£  and  e„  are  used  and  a  simple  sum  formed  to 
account  for  both  effects,  as  suggested  by  Patel  and  Lee. 
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3.  No  lag  equation  is  used  in  conjunction  with  e^  or  e^  because  none 
could  be  found  that  was  universal  for  the  range  of  geometries  tested, 
and  further,  to  separate  out  the  x-dependency  requires  a  separate  lag 
equation  for  e^/u  and  e^/u  with  the  remaining  y  dependency  omitted. 

4.  At  a  particular  x-station,  the  thick  boundary  layer  eddy  viscosity 

(without  extra  rates  of  strain)  is  computed  according  to  Eqs.  (11), 

(12)  and  (24)  and  y ^  found.  Then  the  mixing  length  ratio  l/l  is 

computed  from  Eq .  (37)  for  0  <  y  <  y  .  From  Eq .  (11),  (e  )  ~  £  ^ 

C  O  1  o 

and  hence  the  inner  eddy  viscosity,  modified  for  extra  rates  of  strain, 
is  given  by 


ei  =  A  (eo).  ,  0  <  y  <  yc  , 


(43) 


SL 

where  A  =  y-  is  given  by  Eq.  (37)  and  (e  ) .  is  given  by  Eqs.  (11)  and 
x  o  i 

o 

(12).  For  Yc  <  y  <  6  the  value  of  A  at  y  =  y^  is  used,  denoted  by  Ac« 
Thus  the  modified  outer  eddy  viscosity  is 


e  =  A  (e  )  ,  y  <  y  <  6  , 

o  c  o  o  c  -  - 


(44) 


and  (e  )  is  given  by  Eq.  (24). 
o  o 


Results 

The  boundary  layers  on  three  bodies  of  revolution  each  having  a  different 
longitudinal  curvature  history,  have  been  calculated  with  the  ARL  finite 
difference  boundary-layer  code  using  the  foregoing  algebraic  turbulence 
model  with  and  without  extra  rates  of  strain.  For  each  of  these  bodies  a 
sufficient  amount  of  detailed  experimental  measurements  exists  with  which 
to  make  comparisons.  The  bodies  are: 

1.  F-57  low  dTag  body,  ref.  1 

2.  Modified  spheroid,  ref.  6 

3.  NSRDC  body  (afterbody  no.  1),  ref.  7 
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The  geometries  are  shown  in  Figs.  6-8. 

Boundary-layer  input  parameters  for  the  three  cases  are  listed  in 
Table  1.  Transition  locations  in  the  calculations  were  chosen  so  that  0 
in  the  laminar-turbulent  region  closely  fitted  the  experimental  data. 

For  the  step  size  in  the  direction  normal  to  the  body  surface  a  geometric 
progression  was  used  where  the  transformed  step  size  ratio  was  taken  to  be 
1.10.  At  the  stagnation  point  boundary-layer  infinity  was  taken  to  be 
=  5.0,  where  f)  is  the  Mangler-Levy-Lees  variable  defined  by  Eq .  (7.3.4) 
of  ref.  2.  The  pressure  distribution  in  each  case,  obtained  by  the  Douglas 
Neumann  procedure,  was  modified  in  the  stern  region  to  coincide  with  the 
experimental  measurements.  All  cases  were  run  on  a  VAX  11/780  computer  using 
double  precision  arithmetic. 

The  tuning  of  the  modified  turbulence  model  was  performed  using  the 
F-57  and  modified  spheroid  data  only.  The  best  agreement  with  these  data 
was  found  when  the  constant  a  in  Eq.  (38)  had  the  value  3.0.  The  solution 
for  the  NSRDC  body  was  then  run  using  this  value  of  a.  In  all  three  cases 
solutions  were  also  obtained  using  the  unmodified  turbulence  model  (a  =  0). 

The  reason  for  using  the  F-57  and  modified  spheroid  in  the  turbulence 

model  tuning  process,  aside  from  the  high  quality  of  the  data,  is  that  these 

bodies  have  very  different  longitudinal  curvature  histories.  Consequently, 

the  contributions  of  e  and  e.  to  the  individual  boundary  layers  will  be 

t 

different. 

For  the  F-57  and  modified  spheroid,  mean  velocity  profiles,  Reynolds 
stress  profiles,  skin  friction  coefficient,  boundary  layer  shape  factor  and 
momentum  deficit  area  are  presented,  whereas,  for  the  NSRDC  body  only  mean 
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velocity  profiles  and  skin  friction  coefficient  are  given.  Comparisons  of 
calculations  with  experiment  for  the  F-57  body  appear  in  Figs.  9-17,  for 
the  modified  spheroid  in  Figs.  18-28  and  for  the  NSRDC  body  in  Figs.  29-32. 

For  the  F-57  body  calculations  near  the  stern  using  the  modified 
turbulence  model  are  in  better  agreement  with  experiment  for  the  mean  velocity 
and  Reynolds  stress  profiles  than  corresponding  calculations  with  the 
unmodified  model  (see  Figs.  9-14).  At  the  last  station,  xq  =  0.96,  the 
modified  turbulence  model  suppresses  the  Reynolds  stress  too  much  near  the 
wall  and  produces  an  overcorrection  in  the  velocity  profile.  Nevertheless, 
the  modified  result  is  still  in  better  agreement  with  experiment  than  the 
unmodified  one.  Suprisingly  enough,  the  shape  of  the  Reynolds  stress 
profile  is  qualitatively  correct  for  the  modified  model. 

The  modified  model  reduces  the  skin  friction  coefficient  slightly 
below  the  measured  values,  as  shown  in  Fig.  15,  indicating  that  £/£.q,  as 
given  by  Eq.  (37),  is  decreasing  too  rapidly  near  the  wall.  One  possible 
way  of  suppressing  this  rapid  drop  off  would  be  to  use  a  damping  factor  in 
Eq.  (37),  similar  to  Eq.  (13).  As  seen  in  Figs.  (16)  and  (17),  the  shape 
factor  prediction  is  improved  somewhat  by  the  modification  while  the 
momentum  deficit  area  remains  almost  unchanged,  as  would  be  expected. 

In  the  case  of  the  modified  spheroid,  for  mean  velocity  and  Reynolds 
stress,  as  shown  in  Figs.  19-25,  the  modified  turbulence  model  gives  results 
in  better  agreement  with  experiment  than  those  of  the  unmodified  model. 

Here  the  improvement  is  not  quite  as  dramatic  as  in  the  case  of  the  F-57 
body,  except  at  the  last  station,  =  0.99.  The  agreement  can  be  improved 

by  increasing  a  to  about  3.5  or  4.0,  but  then  the  results  for  the  F-57 
body  become  worse  because  the  correction  becomes  too  large  near  the  tail. 
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In  Figs.  26  and  27  the  modification  is  seen  to  improve  agreement  for 
the  skin  friction  coefficient  and  shape  factor.  The  momentum  deficit  area, 
shown  in  Fig.  28,  is  essentially  unchanged  by  the  modification,  as  was  the 
case  for  the  F-57. 

The  mean  velocity  profile  measurements  given  in  ref.  7  for  the  NSRDC 
body  are  in  cylindrical  coordinates  (xq,  r)  with  corresponding  velocity 
components  (u^,  v  ).  To  obtain  calculated  velocity  profiles  in  cylindrical 
coordinates  required  the  use  of  double  quadratic  interpolation  plus  the 
usual  axis  rotation.  The  mean  velocity  profile  comparisons  are  shown  in 
Figs.  29-31.  As  in  the  previous  two  cases,  the  modified  turbulence  model 
gives  results  which  are  in  very  good  agreement  with  the  experimental  data 
and  are  superior  to  the  results  of  the  unmodified  model.  Finally,  in  Fig.  32 
the  modification  is  seen  to  lower  the  skin  friction  coefficient  as  much  as 
20  percent  below  the  measured  values.  This  situation  is  the  same  as 
occurred  for  the  F-57  body  (see  Fig.  15). 

Conclusions 

1.  The  inclusion  of  extra  rates  of  strain  in  a  two  piece  algebraic 
turbulence  modex  has  been  found  to  provide  a  definite  improvement  in  the 
prediction  of  mean  velocity  and  Reynolds  stress  profiles  as  well  as  shape 
factor  in  the  tail  region  of  an  axisymmetric  boundary  layer. 

2.  Both  longitudinal  and  transverse  extra  rates  of  strain  must  be 
incorporated  in  the  turbulence  model  to  produce  best  agreement  with 
experiment.  In  the  present  treatment  the  equilibrium  values  of  e  and 

e.  are  added  together  to  form  an  effective  e  which  is  then  used  in  a  simple 
nonlinear  length  scale  ratio  formula  to  modify  the  mixing  length.  This 
formulation,  a  modification  of  a  scheme  proposed  by  Patel  and  Lee,  has 
been  found  to  be  the  best  of  several  tried. 
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3.  The  best  value  of  the  empirical  constant  a  in  Eq .  (37)  for  the  length 
scale  ratio  has  been  found  by  numerical  experiments  to  be  3.0. 

4.  The  predicted  skin  friction  coefficient  distribution  in  the  vicinity  of 
the  tail,  as  given  by  the  turbulence  model  modified  for  extra  rates  of 
strain,  has  been  found  to  lie  generally  below  measured  values.  This 
defect  appears  to  be  related  to  the  overly  rapid  decay  of  length  scale 
ratio  near  the  wall.  A  possible  remedy  is  to  employ  a  damping  function, 
similar  to  Eq .  (13),  in  the  vicinity  of  a  wall. 

The  axisymmetric  boundary-layer  computer  code  with  the  modified 

turbulence  model  described  in  this  report  is  available  from  the  author. 
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Figure  7.  Modified  Spheroid  Geometry. 
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Figure  8. 


NSRDC  Body  Geometry. 


Figure  9.  Mean  Velocity  Profile  at  x  =  0.88,  F-57  Body. 


Figure  11.  Mean  Velocity  Profile  at  x  =  0.96,  F-57  Body 


Figure  15.  Skin  Friction  Coefficient  Distribuiton,  F-57  Body, 
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Figure  17.  Momentum  Deficit  Area  Distribution,  F-57  Body 
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